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ABSTRACT 

We present a comparison of CN bandstrength variations in the high-metallicity globular clusters NGC 
6356 and NGC 6528 with those measured in the old open clusters NGC 188, NCG 2158 and NGC 
7789. Star-to-star abundance variations, of which CN differences are a readily observable sign, are 
commonplace in moderate-metallicity halo globular clusters but are unseen in the field or in open 
clusters. We find that the open clusters have narrow, unimodal distributions of CN bandstrength, 
as expected from the literature, while the globular clusters have broad, bimodal distributions of 
CN bandstrength, similar to moderate-metallicity halo globular clusters. This result has interesting 
implications for the various mechanisms proposed to explain the origin of globular cluster abundance 
inhomogeneities, and suggests that the local environment at the epoch of cluster formation plays a 
vital role in regulating intracluster enrichment processes. 

Subject headings: Globular clusters: individual (NGC 6356, NGC 6528) - Open clusters: individual 
(NGC 188, NGC 2158, NGC 7789) - Stars: abundances 



1. INTRODUCTION 

Intracluster abundance variations are well-known 
among globular cluster stars on the red giant branch, 
subgiant branch, and main sequence (see Gratton et al. 
2004 for a thorough review). Within halo globular clus- 
ters, at fixed stellar luminosity, some stars are observed 
to have a typical Population II abundance pattern, while 
others are depleted in carbon, oxygen, and magnesium, 
and enriched in nitrogen, sodium, and aluminum. This 
latter abundance pattern is most easily observed as un- 
usually strong CN absorption in low-resolution spectra, 
whereby it was first discovered. Consequently such stars 
are referred to as "CN-strong," and the entire abundance 
pattern from carbon through aluminum is implied by 
that name. The existence of CN-strong stars in globular 
clusters is presently interpreted as a sign of an early en- 
richment process whose details are murky (e.g.. Cannon 
et al. 1998, Bekki et al. 2007). This process appar- 
ently occurs in all halo globular clusters (|Gratton et al.l 
I2004f) . but not in open clusters (e.g., N orris fc Smith 
1985 , Smith & Norris 1984), the halo field (jCratton et all 
[200nh or in the dwarf sph eroidal satellites of the Milky 
Way (jShetrone et al.ll2003l ). 

There are several possible explanations for the exis- 
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tence of CN-strong stars in globular clusters, includ- 
ing pollution of cluster-forming gas by AGB star winds 
(e.g., Cottrell & Da Costa 1981), an early generation of 
stars with a top-heavy mass function (e.g.. Cannon et al. 
1998), or Wolf-Rayet stars (e.g.. Smith 2006). There is 
also an emerging paradigm that sufficiently massive glob- 
ular clusters form their stars in two generations (e.g., 
Carretta et al. 2008 and D'Ercole et al. 2008). In 
this picture, cluster mass loss is driven by supernovae, 
primarily affecting the first generation of stars and ac- 
counting for the varying ratios of halo-like to enriched 
stars and differing horizontal branch morphologies. All 
of the various models have a few common characteristics: 
the process must happen quickly, since there is a limited 
time from the first star formation in a globular cluster 
until the O.8M0 stars that are presently red giants have 
finished forming. All models rely on intermediate- to 
high-mass stars as a source for polluting material, both 
because of their relatively short lifetimes and because 
the existence of anticorrelated Mg and Al abundances 
implies high-temperature fusion. All models also require 
a proto-cluster to have a significant gravitational poten- 
tial (to retain enriched gas) and a sufficiently rarefied 
environment (to prevent ram pressure stripping). Given 
these characteristics, it can be understood qualitatively 
why old, isolated halo globular clusters would contain 
CN-strong stars while lower-mass disk objects like open 
clusters would not. 

Disk globular clusters, with their relatively high metal- 
licities, provide an interesting intersection of those two 
populations, and a convenient testbed for models of the 
primordial enrichment process. If the process is prohib- 
ited by in-disk orbits or high metallicity, then disk glob- 
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ular clusters will have homogeneous abundance patterns 
shnilar to those in open clusters. However, if the high 
initial mass and relatively deep potential well of a proto- 
globular cluster is the key to allowing enrichment, disk 
and bulge globular clusters could also contain a subpopu- 
lation of stars with the same CN-strong abundance pat- 
tern observed in halo globular clusters (if the yields of 
the enriching stars scale with metallicity) . 

It is in this context that we consider our data set, 
which is comprised of Keck/LRIS blue spectra of 24 up- 
per red giant branch stars in the high-metallicity disk 
globular cluster NGC 6356 and 11 similar stars in NGC 
6528. NGC 6356 is in a location consistent with being 
a disk globular cluster (e.g., Zinn 1985): 2.7 kpc above 
the Galactic plan e, on the far side of the Galactic center 
(|Bica et al.lll994[ ). Its metallicity is [Fe/Hj^ - 0.50 ac- 
cording to the 2003 revision of the lHarrisI (|1996[) catalog. 
NGC 6528 is located in Baade's Window, 0.6 kpc from 
the Galactic center and 0.6 kpc below the plane. As such , 
it is considered a bulge globular cluster (Minniti ' 1995f ). 
with Zoccali et al. (2004) and Origha et al. (2005) find- 
ing metallicities of [Fe/H]= —0.1 and —0.2 dex from opti- 
cal and infrared spectroscopy, respectively. We compare 
these spectra to Lick/Kast blue spectra of 22 red giants in 
the old open clusters NGC 188 ([Fc/H]= -0.06, Carraro 
& Chiosi 1994), NGC 2158 ([Fe/H]= -0.60, Carraro et 
al. 2002), and NGC 7789 ([Fe/H]= -0.04, Tautvaisiene 
et al. 2005). These are all northerly clusters, with esti- 
mated ages ranging from 1.3 Gyr for NGC 7789 to 7.5 
Gyr for NGC 188 (both age estimates from Carraro & 
Chiosi 1994). We measure the strength of the CN ab- 
sor ption band at 4215A with the index 5(4142) defined 
bv iNorris fc FreemanI (|1979l ). This quantity, like most 
spectral indices, measures the magnitude difference be- 
tween the integrated flux inside the absorption feature in 
question and in a nearby continuum region. By consid- 
ering CN bandstrength variation as a proxy for the full 
enriched abundance pattern, we hope to explore whether 
NGC 6356 and NGC 6528 exhibit CN bimodality like 
their low-metallicity halo counterparts do, or whether 
their similarity to open clusters includes a homogeneous 
abundance distribution along with near-Solar metallici- 
ties and disk orbits. 

2. THE DATA SET 

Spectra in NGC 6356 were acquired on two nights in 
2002 June, using the Low Resolution Imaging Spectrom- 
eter (LRIS, Oke et al. 1995) at the Keck I Telescope on 
Mauna Kea. Spectra in NGC 6528 were taken on two 
nights in 2003 June with the same instrumental setup. 
Using a mirror in place of a dichroic, all light was directed 
to the mosaic of two 2Kx4K Marconi CCDs in the blue 
camera. The 400/3400 grism and a slit width of 1" (8.7" 
for flux standard stars) resulted in a pixel spacing of 1.1 A 
and a resolution of 8A. Typical exposure times of 1200- 
1500 seconds produced a signal-to-noise ratio (S/N) per 
pixel of roughly 60 at 4320A, just redward of the G band. 

The open cluster stars were observed with the Shane 
3m telescope at Lick Observatory on two nights in 2005 
February, three nights in 2005 September, and three 
nights in 2006 August. In each case a mirror was used to 
direct all light to the blue side of the Kast double spec- 
trograph. The 600/4310 grism and a slit width between 
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NGC 6356: Star identifiers and photometry 
are taken from ISandaee fc Wallersteiiil 1 IT9BQ ).'' NGC 
6528: Star identifiers and photometry are from 
Ivan den Bereh fc Youneeri l lT979h . 

1" and 2", depending on weather (9" for flux standard 
stars), produced a pixel spacing of 1.8A and a resolution 
of 5.4A. Exposure times varied widely because of a range 
of two magnitudes in distance modulus, and the typical 
signal-to-noise ratio at 4320A is 150. Tabled lists iden- 
tifying information, photometry, exposure time and the 
CN bandstrength index S'(4142) for each globular cluster 
giant observed. Analogous information for the open clus- 
ter giants, and membership probabilities where available, 
are given in Table [21 

Data reduction was accomplished for both the Keck 
and Lick data sets using the XIDL "low-redux" pro- 
grams^. The instrument-specific programs handle bias 
subtraction, flat-fielding, cosmic ray removal, object 
identification and extraction, sky subtraction, flexure 
correction, wavelength calibration, atmospheric correc- 
tion, coadding, and flux calibration. The standard stars 
Feige 34, Feige 110, BD +28 4211, and BD -^33 2642 were 
used for flux calibration, with "tr ue" sp ectra taken from 
the CALSPEC database (Bohhi]|l99i). Figure [D shows 
a typical flux-calibrated spectrum from each of the Kast 
(open cluster) and LRIS (globular cluster) data sets, with 

^ Produced by J.X. Prochaska and available at 
|http://www. ucoIick.org/~xavier/LowRedux/ 1 
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^ NGC 188: Star identifiers and photometry are taken from ISandage| (11963), 
proper motion-based membership probabilities are fromlStetson ct al. (2003), 
radial velocity-based membership probabilities are from [Gcllcr ct al. (2003)- 
^ NG C 2158: star identifiers and photometry are taken from lArp fc CuffevI 

1196^. " NGC 778 9: Star identifiers and photometry are taken from 

IMcNamar a fc Solomon f 1981), radial velocity- based membership probabilities 
are from[Gim ct al. (1998). 



the science band (4120AA to 4216AA) and comparison 
band (4216AA to 4290AA) of S'(4142) shown as regions 
with closely spaced shading and widely spaced shading 
respectively. 

2.1. Cluster Membership 

Because all of the star clusters observed for this study 
orbit within the disk or bulge, the line of sight to each 
of them contains more foreground stars than the line of 
sight to a typical halo globular cluster. Because of the 
near-solar metallicities of the cluster stars and the typ- 
ically low systemic radial velocities of the clusters rela- 
tive to the Sun, two of the most straightforward meth- 
ods for identifying foreground stars in low-metallicity 
globular cluster fields are ineffective with the present 
data set. Three of our five target lists (for NGC 6356, 
NGC 6528, and NGC 2158) were constructed from photo- 
graphic studies without the astrometric or spectroscopic 
information necessary to identify foreground late-K and 
M dwarfs. 

We use a combination of visual spectrum inspection, 
photometry, and spectral indices focused on gravity- 
sensitive features to identify foreground dwarfs and other 
non-useful stars in our data set and remove them from 
further analysis. We reject NGC 6356 17 because of TiO 
bands not seen in any of the other NGC 6356 stars, and 
NGC 6356 51 because it is signi ficantly bluer and has 
a distinctly smaller Ca K index (jMorrison 
than all other stars observed in NGC 6356. We also 
drop NGC 6528 11-27 from the sample beca use its A3883 
CN b and, measured by the index S'(3839) (iBrilev et al.l 
|199C1| ). is abnormally small. In NGC 188, though the stars 
are all confirmed as cluster members by pr oper motion 
(jStetson et al.l |2004[ ) and radial velocity (|Geller et al.l 



l2008f ). star 1 001 is thought to be a po st-coalescence FK 
Comae star (Harri s fc McClurd fl985[ l and star 3018 is 
significantly brighter than the rest of the sample. We are 
therefore unsure whether 3018 is an AGB star, and we 
remove it from the analysis. NGC 7789 also has proper 
motion-based me mbership probabilities in the literature 
(|Gim et al.lll998f ). We do not see any outliers in the pho- 
tometry or measured indices in NGC 7789, and include 
aU of the NGC 7789 stars hsted in Table [2] in the analy- 
sis. In all, we exclude 5 of the 62 stars initially observed 
from the final sample, leaving 24 stars in NGC 6356, 11 
in NGC 6528, 6 in NGC 188, 6 in NGC 2158, and 10 in 
NGC 7789. 

3. CN BANDSTRENGTH DISTRIBUTION 

We measure the CN bandstrength index S'(4142) from 
each of our fiux-calibrated spectra. Because it is de- 
fined as the magnitude difference between the integrated 
flux inside the 4215A CN band and a nearby CN-free re- 
gion, larger values of 5(4142) reflect stronger CN absorp- 
tion. We determine the error on the index measurement 
on a cluster-by-cluster basis. We define a pseudoindex 
5*^(4142), with the same comparison band as the original 
S'(4142) and the science band shifted to 4027A through 
4077A, a CN-free region of the spectrum. The scatter in 
S'p(4142) at a given absolute magnitude will then be inde- 
pendent of CN bandstrength, and will be representative 
of the scale of the random errors in S'(4142). Figure [2] 
shows 5(4142) (upper panel) and 5p(4142) (lower panel) 
versus My for our NGC 7789 stars. As expected from 
previous abundance-inhomogeneity studies of open clus- 
ters, the stars stay within a fairly small range in 5(4142). 
In NGC 7789 the scatter in 5p(4142), which is 0.05, is 
comparable to that in 5(4142), implying that there is 



4 




Fig. 1. — Upper panel: representative spectra from the data set. 
The upper spectrum is NGC 7789 star 614, taken with the Kast 
spectrograph at Lick Observatory. The lower spectrum is NGC 
6356 111, taken with the LRIS spectrograph at Keck Observatory. 
The closely spaced shaded region shows the science band of the 
index 5(4142), and the more widely spaced shading marks the 
5(4142) comparison band. Lower panel: comparison of CN-strong 
and CN-eak spectra in the region of the 4215A CN band. The 
thinner line is the spectrum of NGC 6356 113, which has a stronger 
CN band than the heavier line, which is the spectrum of NGC 6356 
111. 

very little CN abundance variation among the stars. It 
is also worth noting that 5*^(4142) rises with brightening 
My in Figure 2, while S'(4142) does not. This could be 
a result of the larger wavelength range spanned by the 
bandpasses of 5'j,(4142), which makes it more sensitive 
to changes in stellar temperature. 

This pseudoindex error, calculated as the RMS of 
5*^(4142) about a best-fit fine in the S'p(4142) - My plane, 
is calculated for each of the clusters in our sample indi- 
vidually. For NGC 6356 the scatter in 5^(4142) is 0.048 
and for NGC 6528 the scatter is 0.082; a representative 
value of ±0.05 magnitudes is adopted as the error in all 
measured values of S'(4142). After inspecting some of 
the spectra, we feel that large-scale variations in con- 
tinuum shape are most likely to be responsible for the 
variations in S'p(4142). These may be a result of at- 
mospheric diffraction causing the scattering of blue light 
preferentially out of the slit, particularly in NGC 6356 
and NGC 6528, which are fairly southerly for observa- 
tion from the Keck I telescope. Since the bandpasses of 
S'p(4142) are more widely separated in wavelength than 
those of S'(4142), and extend further into the violet, the 
former index will be more sensitive to observational ef- 
fects such as differential light loss at the spectrograph 
slit, and errors in flux calibration. As such, our approach 
may overestimate the observational error in 5(4142) for 
the globular cluster observations. 

The left panel of Figure [3] shows the CN bandstrength 
index 5(4142) versus absolute visual magnitude for NGC 
6356. There is clearly a vertical range in 5(4142) at 
each My, which is also seen in every other CN-bimodal 
globular clu s ter. A dapting the analysis method of the 
iNorris et af] ()1981f ) study of CN variations in NGC 6752, 
we fit a line to the lower collection of points and mea- 
sure the quantity (55(4142), the vertical distance from 
each point to the line. In this study we set the slope 



Fig. 2. — Values of the CN index 5(4142) (upper panel) and 
the pseudoindex 5p(4142) (lower panel) measured for red giants in 
NGC 7789. As described in the text, the RMS scatter in 5p(4142) 
is used for estimation of measurement error on the index 5(4142). 

of the baseline to zero, in keeping with previous studies 
on 4215A CN bandstrength in high-metallicity globular 
clusters (e.g., Briley 1997). The right panel of Figure [3] 
shows a generalized histogram of (55(4142) values mea- 
sured from the left panel. This curve is produced by rep- 
resenting each (55(4142) point as a Gaussian with a width 
(T equal to the measurement error, which in this case is 
0.05. The sum of those individual Gaussians is then a re- 
alistically smoothed histogram. Figure S] shows the anal- 
ogous data in NGC 6528, with a generalized histogram 
calculated in the same way. With fewer data points and 
a smaller gap in 5(4142) between relatively CN-weaker 
and CN-stronger groups, the generalized histogram for 
NGC 6528 looks less bimodal than the generalized his- 
togram for NGC 6356. The dashed curve in the right 
panel of Figure [4] shows the generalized histogram that 
results if a is set to 0.08, the measured standard devia- 
tion in 5p(4142) in NGC 6528. Increasing a smooths the 
curve dramatically, reducing what appears in the solid 
curve to be a distinct CN-stronger group to a broad, 
asymmetric distribution. 

Figure [5] shows generalized histograms of (55(4142) for 
the open clusters NGC 188, NGC 2158, and NGC 7789, 
all calculated in the same way as for the globular clusters, 
with flat baselines and a set to 0.05. The curves are all 
normalized at the peak for easier intercomparability, and 
the generalized histogram for NGC 6356 is overplotted in 
each panel as a dotted curve. It is clear that the distri- 
bution of CN bandstrength in the three open clusters is 
single-peaked, and that there are two groups in CN band- 
strength in NGC 6356, but NGC 6528 is more difficult 
to categorize. It is possible that the difference in [N/Fe] 
between CN-enhanced and CN-normal stars is smaller in 
Solar-metallicity globular clusters, since the total mass 
in nitrogen required to produce an abundance gap of up 
to 1 dex is far larger than in low-metallicity globular 
clusters. The broad but not clearly bimodal distribution 
of CN bandstrengths in NGC 6528 could be a sign that 
there is an upper limit on the amount of nitrogen that 
can be injected into feedback material in high-metallicity 
globular clusters. 
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Fig. 3.— Measured values of the CN index 5(4142) in NGC 
6356 {lejt panel) and a generalized histogram of the offset <5S'(4142) 
between the data points and the horizontal baseline {right panel). 



4. RESULTS 

The homogeneity of CN bandstreng ths in the open 
cluste rs surveyed is an expecte d resu lt: iNorris fc Smith 
' 1985|) . iSmith fc NoTrii (|1984f ). and iTautvaisiene et al 
2005f ) have all demonstrated small ranges in band- 



strengths or abundances in these clusters. However, the 
disk/bulge globular clusters NGC 6356 and NGC 6528 
have distributions of CN bandstrength that clearly re- 
semble those observed in halo globular clusters like NGC 
6752 and M3, and that is surprising. We interpret the 
presence of CN-strong stars in NGC 6356 and NGC 
6528 to mean that the same primordial enrichment pro- 
cess took place in disk and bulge globular clusters as in 
halo globular clusters. Higher resolution spectroscopy is 
needed to determine whether the full "enhanced" abun- 
dance pattern from carbon through alumin um is present 
in the CN-strong stars. On a related note. [Origlia et all 
([200^ have found evidence of a Mg-Al abundance anti- 
correlation, such as often accompanies CN variations in 
globular clusters, in the bulge globular NGC 6441 ([Fe/H] 
= -0.5). 

The presence of CN-strong stars in disk and bulge 
globular clusters provides strong constraints on the pri- 
mordial enrichment process. Since the Galactic orbits of 
NGC 6356 and NGC 6528 do not allow for long periods 
of time outside the disk or bulge, the transfer of enriched 
material from evolved higher-mass stars to still-forming 
low-mass stars must have been both relatively fast and 
efficient in order to avoid ram pressure stripping from 
a primordial environment that was likely much denser 
than at present. In addition, the absolute amount of ni- 
trogen needed to create a difference of between 0.3 and 
1.0 dex in [N/Fe] between CN-weak and CN-strong stars 
is significantly larger at near-Solar metallicity than it is 
at typical halo metallicity. This implies that the material 
transferred to low-mass stars must have been more en- 
riched in nitrogen in the higher-metallicity globular clus- 
ters. Such a requirement may pose a problem for AGB 
star enrichment scenarios. The 4 and 5 Mp^ AGB mod- 
els wit h Reimers mass loss rates of iKarakas fc Lattanziol 
(|2007f) eject a total mass of ^^^N of 7.7 x lO'^ Mq and 



Fig. 4. — Measured values of the CN index S(4142) in NGC 
6528 {lejt panel) and a generalized histogram of the offset <5S'(4142) 
between the data points and the horizontal baseline {right panel). 
The presence of a smaller gap between relatively CN-stronger and 
CN-weaker groups than in NGC 6356, and the smaller number of 
CN-stronger stars, create a less bimodal-looking generalized his- 
togram than in Figure |3] The dashed curve shows the generalized 
histogram calculated with a cr of 0.08. 



3.6 X 10-2 Mq respectively for Z = 0.0001, 2.3 x lO'^ 
Mq and 6.3 x 10"^ M© respectively for Z = 0.004, 
and 2.2 x 10"^ Af© and 5.4 x 10"^ M© respectively for 
Z = 0.008. There is no indication from these models 
that early AGB stars in bulge or disk globular clusters 
would eject a proportionally larger amount of nitrogen 
than intermediate-mass AGB stars in halo clusters. To 
obtain comparable enhancements of [N/Fe] in CN-strong 
giants in both disk and halo globular clusters, it may be 
necessary to assume that the former were enriched by a 
larger number of AGB stars per unit cluster mass than 
halo clusters. 

The evidence in this article indicates that disk and 
bulge globular clusters contain a population of CN- 
enhanced giants, presumably analogous to those in halo 
globular clusters, whereas open clusters of similar metal- 
licity do not. This may indicate that open clusters were 
not massive enough to experience early enrichment of 
the type evinced by globular clusters. Alternatively, disk 
and bulge globular clusters may have had a larger than 
typical number of AGB stars per unit mass. Could this 
have been a consequence of the environment in which 
the bulge and disk globulars formed? Conditions in 
the central parts of the proto-Galaxy may have favored 
high- and intermediate-mass star formation under star- 
burst conditions. Open clusters, forming further out in 
the Galactic disk in more quiescent molecular clouds, 
could have resulted from star formation activity more 
conducive to lower-mass stars. 

Indeed, this consideration of AGB-star yields suggests 
that self-enrichment of disk and bulge globulars may been 
have associated with more top-heavy stellar mass func- 
tions than were present even in halo globular clusters. 
Many halo globular clusters may have formed within 
proto-dw arf galaxies or massiv e clouds of the type sug- 
gested by Searle fc ZiniJ ()1978l ). and their chemical evo- 
lution would have been set by the environment within 
such systems. Recent studies of multiple stellar popu- 
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Fig. 5. — Generalized histograms of the ON bandstrength index 
(55(4142) for each open cluster in this survey, with the generalized 
histogram of <5S'(4142) in NGC 6356 overplotted as a dotted line. 
The open clusters all have single-peaked distributions in CN band- 
strength, in clear contrast to the disk globular cluster NGC 6356. 



lations i n Galactic and extr a galactic globular clu sters, 
such as iMilone et all (|2008f) . iPiotto et all (|2007f) . and 
iMackevet all (|2008[ ). make it clear that their star forma- 
tion took place in a complex and dynamic environment. 
The sites of disk and bulge cluster formation may have 
been strongly influenced by conditions within the chaotic 
and violent inner regions of the proto-Galaxy. Cloud col- 
lisions may have been fairly common in the still-forming 
bulge, resulting in starbursts that were relatively efficient 
in high- and intermediate-mass star formation, leading to 
pronounced chemical enrichment that built up not only 
high overall metallicities, but also high [N-Na-Al/Fe] ra- 
tios as well. As such, the study of abundance inhomo- 
geneities in disk and bulge globular clusters could offer 
the promise of probing some of the most efficient primor- 
dial environments for chemical evolution. 
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